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We revisit recent developments in the theoretical foundations of time-dependent density functional 
theory (TDDFT). TDDFT is then applied to the calculation of total cross sections for ionization 
processes in the p-He collision system. The Kohn-Sham potential is approximated as the sum of the 
Hartree-exchange potential and a correlation potential that was proposed in the context of laser- 
induced ionization. Furthermore, some approaches to the problem of calculating the ionization 
probabilities from the density are discussed. Small projectile energies (<5keV) are considered as 
well as those in the range from 5 to 1000 keV. Results are compared with former calculations and 
with experimental data. We find that the correlation potential yields no obvious improvement 
of the results over the exchange-only approximation where the correlation potential is neglected. 
Furthermore, we find the problem of calculating the desired observables crucial, introducing errors 
of at least the same order of magnitude as the correlation potential. For the case of small energies 
we find that trajectory effects play an important role: the ionization cross sections are enlarged 
significantly if curved instead of straight-line trajectories are used for the projectile motion. 
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I. INTRODUCTION 

There are several good reasons for studying atomic col- 
lisions. Probably the most important one on the funda- 
mental side is their suitability to shed light on the quan- 
tum dynamical many-particle problem. Accordingly, 
they often serve as test-beds for many-particle methods. 
This is especially true for the antiproton-helium system 
under consideration in the present paper, which is ar- 
guably the simplest collision system with more than one 
electron (Complications associated with electron trans- 
fer processes or exchange effects are absent. This is in 
contrast to bare ion and electron impact.) This makes 
it a problem in which on the one hand the interaction 
between the electrons has to be taken into account, but 
which, on the other hand, is still simple enough to enable 
a numerical solution of the time-dependent Schrodinger 
equation with present-day computers. These solutions 
can then serve as benchmarks for other methods, e.g., 
time-dependent density functional theory (TDDFT). 

In the present paper we explore the TDDFT approach 
to the problem at hand and are especially concerned with 
correlation effects. There are two ways in which they are 
relevant in TDDFT: First, an approximation for the cor- 
relation potential is necessary to calculate the time evo- 
lution of the density (called dynamical correlation in the 
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following). Second, the observables of interest must be 
expressed as density functionals (we will call this func- 
tional correlation in the following). We investigate a re- 
cent model suggested in the context of laser-induced ion- 
ization [1] for the treatment of the dynamical correlation, 
and one [2] for the functional correlation, and we give an 
estimate for the error induced by functional correlation. 
Furthermore, we will apply these approximations to col- 
lisions at lower energies than are usually considered. 

In Sec. II, after giving a short introduction to the fun- 
damentals of TDDFT, we introduce the time-dependent 
exchange-only (hereafter: x-only) approximation investi- 
gated in Ref. [3] and the above-mentioned approxima- 
tions for correlation effects. A few remarks on the com- 
putational methods used are also provided. In Sec. Ill we 
present our results and compare them with experimental 
data and several other calculations. Conclusions are of- 
fered in Sec. IV. Atomic units are used unless indicated 
otherwise. 



II. THEORY AND COMPUTATIONS 

A. Time-dependent density functional theory 
(TDDFT) 



TDDFT is a method for the investigation of systems of 
N indistinguishable, interacting particles, described by a 
time-dependent Schrodinger equation (TDSE) 

i$|*(t)> = #(*)!*(*)> V t g [0,7], |tf (0)> = |* > (1) 
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with the Hamiltonian 



H(t) = T + W + V(t), 



(2) 



where T is the kinetic energy operator, W is the interac- 
tion potential and V an external potential that depends 
on time. We denote by r\ the position, by o\ the spin of 
the Z-th particle and by x\ the combination of both, and 
we set x = (xi . . . Xn). 

The first Runge- Gross theorem [4] states, that if the 
initial state |^o) and the interaction potential W are 
fixed, the time-dependent (one-particle) density function 
n:R 4 ^R+, defined by 

n( ri ,t) = Nj2 J d 3 r 2 ...J d 3 r N \*(x)\ 2 (3) 

determines V up to a time-dependent constant. Hence, 
all physical information is determined as it is invariant 
under the transformation V(t) — >• V(t) + c(t). Especially, 
the state of the system is determined for all times up to a 
merely time-dependent phase. We denote the state that 
yields the density n for the interaction W by |^[W", n]). 

n can be obtained by solving the so-called Kohn-Sham 
equations, which are one-particle TDSEs for noninteract- 
ing particles in an effective, density-dependent, external 
potential V KS [n\: 

idt\tpi(t)) = H KS [n](t)\^(t)), M0)> = |^> (4) 

with H[n] KS (t) = f + V KS [n](t) and the one-particle 
states \(fi) . . . \ipn)- 

Originally (in [4]) V" KS [n] had been derived from the 
Frenkel-stationarity principle as 



V KS [n]( ri ,t) 



V(r u t) + 



6A [W,n] _ SAo[0,n] 
5n(ri,t) Sn(ri,t) 



(5) 



where 



A [W, 



f 

Jo 



dt(V[W,n]\idt-f - W\V[W,n]) (6) 



is the quantum mechanical action. Note that it has re- 
mained a controversial issue if the potential V KS [n](ri, t) 
at a certain time t can depend on the density n(ri, t') at a 
later time t' > t. If that is the case, the problem has to be 
solved by a self-consistency iteration. Recent discussions 
[5] , [6] , [7] were concerned with the question whether such 
an iteration is indeed necessary and whether its conver- 
gence can be guaranteed. In the approximations we use, 
however, this is not an issue as will be seen in Sec. II C. 

As it turned out [8], the definition (5) leads to a vio- 
lation of causality. This means that the derivation from 
the Frenkel-stationarity principle must be incorrect. In- 
terestingly, while the violation of causality has been well- 
known for some time, until recently there was no paper 
that pinpointed the error in the derivation and addressed 
the question whether it can be corrected. This void has 



now been filled by Vignale who showed that the Frenkel 
principle is (in its usual form) not applicable in TDDFT 
[9]. We would like to rephrase the reason why the orig- 
inal derivation of Ref. [4] went wrong: Let be any 
state and let be the solution of the TDSE with the 
Hamiltonian H and a given initial condition. Frenkel's 
stationarity principle states that = if and only if 

• |§(0)> = |*(0)> and |*(T)) = |*(T)) and 

• the action is stationary at = \$). 

However, when one translates this principle from states 
to densities, the validity of the final condition can never 
be ensured: In TDDFT is not known, but only n. 
While the Runge-Gross theorem states that n determines 
|^) it is not true, that n(.,T) determines |\£(T)), instead 
n would have to be known at all times. Hence, the bound- 
ary condition translates to "n(.,t) is the true density of 
the system at all times", making the principle a mean- 
ingless tautology. One is inclined to think that the same 
argument might hold for the initial condition. This is 
true in general, however, if it is additionally assumed 
that |^(0)) is the ground state of the initial Hamilto- 
nian, then ground-state DFT guarantees that n(.,0) de- 
termines |^(0)). 

This explains why the original result of Runge and 
Gross violates causality. It also shows a way out: Leave 
away the final condition. This is exactly what was done 
by Vignale and which leads to the appearance of addi- 
tional terms in the Kohn-Sham potential compared to 
Eq. (5) (see Ref. [9] for details). 



B. The antiproton-helium system 

The system under consideration in the present paper 
consists of an antiproton, passing on a classical path by 
a para-helium atom (where electrons have antiparallel 
spins). The impact parameter is 6, the antiproton's ki- 
netic energy in the laboratory system (the system where 
the helium atom is at rest at the initial time) is E. For 
sufficiently high energies a straight-line trajectory for the 
antiproton can be assumed. At the initial time the (spa- 
tial) Kohn-Sham state can be written as |<E>) = \(f)\(p) (us- 
ing the abbreviation |3>[n]) = |\£[0, n])); the total state is 
antisymmetrical due to the spin function. We neglect all 
spin-dependent interactions, so the spin function can be 
ignored. Since the Kohn-Sham Hamiltonian contains no 
interaction term, the two Kohn-Sham equations of the 
system as well as the corresponding initial conditons are 
identical and the problem is reduced to one and only one 
one-particle TDSE 



H KS [n]Mt)) 



idMt)). 



(J) 



Unfortunately, no exact formula for the Kohn-Sham po- 
tential is known, therefore, approximations are necessary. 
First, the potential is split up: 

V KS [n} = V + V Hx [n} + V c [n}, (8) 
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where V denotes the external potential by the antiproton 
und the helium nucleus. The so-called Hartree-exchange 
potential is defined by 

V Hx [n]{rut) = ^—^ , (9) 

which (in our case) gives 

V Hx [n](r u t) = ^J d 3 r 2 n(r 2 )W(r u r 2 ). (10) 

W is the Coulomb interaction between two electrons, 
i.e., W(ri,r2) = \ Tl -r 2 \ ' ^ e corre l a ^i° n potential V c is 
unknown, and approximations are discussed in the sub- 
sequent section. 
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FIG. 1: Switch function c [Eq. (12)] as a function of the ion- 
ization /. 



C. The correlation potential 

One possible approximation for the Kohn-Sham poten- 
tial (8) is the so-called x-only approximation, where the 
correlation potential is neglected: V c = 0. This was inves- 
tigated in Ref. [3] and will serve as one of our references. 
An explicit model for V c was proposed in the context of 
laser-driven ionization by Lein and Kummel[l]: 



v c (t)= [c(^m)-i]-v Hx {t) 



Here, c is a switch function (see also Fig.l): 

x 



c(x) 



1 + e 



50cc-100 ' 



(ii) 



(12) 



and / is the ionization (average number of electrons in 
scattering states), i.e., 



1 = 2 



lim lim 

R—too t— ^oo 



S(R) 



where S(R) is a sphere around the target with radius 
R. We will call this approximation for the correlation 
potential the LK (Lein and Kummel) approximation and 
the resulting Kohn-Sham potential the LK KS potential. 

The approximation assumes no static correlation (i.e., 
the initial state is identical to the initial state in the x- 
only approximation) and mimics a discontinuity in the 
effective electron-electron potential: when / reaches 1, c 
becomes zero very quickly and therefore 



V Hx + V c = 0. 



(13) 



That means that the electron-electron interaction is 
switched off as soon as (on average) one electron is ion- 
ized. 

Since Vh x is local in time both the x-only and the LK KS 
potentials are local in time and an iteration procedure is 
not necessary. Instead, the Kohn-Sham equation can be 
solved by standard time propagation. 



D. The observable problem 

Even if the correlation potential (and, therefore, the 
whole Kohn-Sham Hamiltonian) were known exactly, the 
calculation of observables would still be a non-trivial 
problem. While the Runge-Gross theorem guarantees 
that the density determines all observables, exact formu- 
las for their calculation are known only in very few cases 
[10]. The observables we are mainly interested in in this 
work are the g-fold ionization probabilities p q+ (b) (depen- 
dent on the impact parameter b) and the corresponding 
total cross sections (TCS) a*+ = / °° dbbp q +(b). They 
fall into the category of not exactly known observables 
(i.e., observables whose functional dependence on the 
density is not known). 



1. Binomial approximation 

Assume that the Kohn-Sham state |<E>) = \(p)\(f) were 
the true state of the system and let p s be the ioniza- 
tion probability of the particle described by \cp). Then 
the two-particle probabilities follow simply from binomial 
statistics as 



P. 



p]+ = 2p a (l-p a ), (14) 

2+ 2 
Ps = Ps' 

As we will see, this a rather rough approximation. 



2. Simplified adiabatic approximation 
One can define the correlation integral 
3 C = lim lim /d 6 r ( p(n, r 2j t) - in(ri, t)n(r 2 , t) ) , 

R—toc t^oo J \ * J 

(15) 
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where p is the two-particle density. With 3 C the exact 
probabilities can be expressed as [2] 



9 1 + 
2+ 



» 0+ + ±a 



(16) 



Since the functional dependence of p on n is unknown, 
3 C cannot be calculated exactly. An approximation was 
introduced by Wilken and Bauer (WB) in Ref. [2]. They 
rewrote 3 C as 

3 C = lim lim d 6 r g c (r u r 2 ,t)n(r u t)n(r 2 ,t) (17) 

R^oo t^oo J 

S(R) 2 

with the correlation function 

p(r u r 2 ,t) 1 



9c(ri,r 2 ,t) 



n(r u t)n(r 2 ,t) 2' 



(18) 



and then replaced n and p in g c by the adiabatic approx- 
imations 



Im + (1 - i> 2 for < I < 1 




for 1 < / < 2 ' 

for < / < 1 
for 1< I < 2 



where ri\ is the ground-state density of He 1+ , n 2 that of 
He and p 2 ist the two-particle ground-state density of He 
[23]. 

We will use a strongly simplified version of this ap- 
proximation: The right panel of Fig. 3 in Ref. [2] shows 
the dependence of 3 C on the ionization / of a laser-driven 
model helium atom for two different sets of laser param- 
eters. In both cases, the dependence is almost identical 
and a polynomial regression yields (see also Fig. 2): 

3 C (I) = -0.352/ 2 - 0.0296/ 3 - 0.0502/ 4 for < / < 1. 

(19) 

In our context the case / > 1 is irrelevant since such 
large ionizations do not arise. We will call the model 
that uses the analytical expression (19) in (16) the sWB- 
approximation. 



3. Average ionization and bounds for the probabilities 

The observable problem can be circumvented by not 
trying to calculate the probabilities, but only the average 
ionization / which can be written as 



n2 + 



(20) 



From Eq. (16) it follows immediately that the ionization 
can be calculated exactly from the non-interacting prob- 
abilities: 



If one wants to benchmark a certain approximation 
to the Kohn-Sham potential, the average ionization is 
the observable that should be used because it can be 
calculated exactly from the density. Furthermore, / is 
obviously an upper bound for p 1+ and for p 2+ (however, 
in our system p 1+ ^> p 2+ , so it is a sensible bound only 
for p 1+ ). There is no known exact lower bound for 
but one can find arguments that for not too high impact 
energies p 1+ > pj + : 



• The WB-approximation as well as the simplified 
version yield 3 C (I) < 0. 

• Eq. (14) implies that < |, but this is not true 
for 

Applying these bounds, the probability can be estimated 
as 



p 1+ = (pI + +pI + )±p 2 s + . 



(22) 



E. Computational details 

The initial condition is the x-only Kohn-Sham ground 
state of para- helium [11]. As mentioned above, the time- 
locality of our approximations of the Kohn-Sham poten- 
tial allows us to solve the Kohn-Sham equation by incre- 
mentally propagating this state. To this end the state 
and the density are expressed in a basis generated by 
the basis generator method [12]. This basis consists of 
bound states (atomic orbitals) up to the principal quan- 
tum number n = 4 and of pseudo scattering states con- 
structed from the bound states as described in Ref. [3]. 
We use the same basis as was used in that work and also 
the same technique and parameters to calculate the av- 
erage ionization / and the Kohn Sham potential in every 
time step. 




I = pi+ + 2p?+. 



(21) 



FIG. 2: Simplified version of the adiabatically approximated 
correlation integral [according to Eq. (19)] as a function of 
the ionization /. 



5 



III. RESULTS 

A. Medium and high energies 

First, we compare our results with experiments by 
Knudsen et al. [13] and Hvelplund et al. [14] (Figs. 3 and 
4). The figures show <j 1+ as a function of the projectile 
energy. In Fig. 3, the observable problem is considered 
by using the bounds according to Eq. (22) and interpret- 
ing pl+ as the best approximation for In Fig. 
4 binomial and sWB approximations are used. 
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FIG. 3: (Color online) Total cross section for one- fold ioniza- 
tion as a function of projectile energy. Present results within 
x-only and LK approximations with error bars denoting the 
bounds according to Eq. (22). Experiments: Knudsen [13], 
Hvelplund [14]. 
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FIG. 4: (Color online) Total cross section for one- fold ioniza- 
tion as a function of projectile energy. Present results with 
various models for the Kohn-Sham potential (x-only, LK) and 
various approaches for the calculation of probabilities (bino- 
mial, sWB). Experiments: Knudsen [13], Hvelplund [14]. 

The x-only approximation and the LK-approximation 
coincide for high energies. This is because for high en- 
ergies there is not enough time for any change in the 
effective electron-electron-potential to make an impor- 
tant contribution. Furthermore, both models coincide 
with the experimental data. For small energies we see 
differences. The correlation enlarges the cross section 



for a simple reason: In the antiproton-helium system 
an ionization larger than 1 is almost never reached, and 
for / < 1 the correlation potential simply amplifies the 
Hartree-exchange potential, which is repulsive. For a 
closer investigation of the correlation potential the he- 
lium atom has to be perturbed more strongly, e.g., by a 
multiply charged ion. This would lead to higher average 
ionizations. 

The effect of the correlation potential on the cross sec- 
tions is of the same order of magnitude as the possible 
deviations induced by the observable problem. Except 
in the range from 10 to 30keV, our results are not in 
contradiction with the results by Knudsen et al. [13] or 
Hvelplund et al. [14]. However, the uncertainties in our 
results due to the observable problem are quite large, so a 
strong agreement cannot be claimed either, which makes 
a more definite statement impossible. 

Supplemental to the comparison with experimental 
data, we compare our calculated TCS to various theo- 
retical investigations of the system in question in Fig. 
5. From the large body of published results we have in- 
cluded only those, which were obtained from a solution of 
the two-electron TDSE, i.e., which include correlation ef- 
fects. A more comprehensive list of previous calculations 
can be found in Ref. [13]. At high energies all calcula- 
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FIG. 5: (Color online) Total cross section for one-fold ioniza- 
tion as a function of projectile energy. Present results within 
x-only and LK approximations with error bars denoting the 
bounds according to Eq. (22). Two-electron calculations: 
Diaz [15], Bent [16], Igarashi [17], Lee [18], Foster [19], Bronk 
[20]. 

tions, except that by Foster et al. [19], coincide. This 
agreement can be explained by perturbation theory. For 
small energies, the calculation by Foster et al. and that 
by Lee et al. agree with our x-only approximation, while 
the one by Igarashi et al. agrees with our LK approxi- 
mation. This leaves the question open, whether the LK 
approximation improves upon the x-only approximation. 

As mentioned in Sec. II D, the average ionization is 
actually better suited for comparisons with experimental 
data and many-particle calculations since it can be ex- 
pressed as a simple density functional and is only sensi- 
tive to the approximation used for the KS potential. Fig- 
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FIG. 6: (Color online) Total cross section for average ioniza- 
tion as a function of projectile energy. Present results within 
x-only and LK approximations. Two-electron calculations: 
Bent [16], Foster [19], Bronk [20]. Experiments: Knudsen 
[13] and [21], Hvelplund [14]. 
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FIG. 7: (Color online) Total cross section for two-fold ioniza- 
tion as a function of projectile energy. Present results with 
various models for the Kohn-Sham potential (x-only, LK) and 
various approaches for the calculation of probabilities (bino- 
mial, sWB). Two-electron calculations: Bronk [20], Foster 
[19], Diaz [15]. Experiments: Knudsen [21], Hvelplund [14]. 



B. Low energies 



ure 6 shows a = <r 1+ + 2<j 2+ as a function of the impact 
energy. Again, good agreement between our calculations, 
experiment and the calculation by Diaz et al. [15] can be 
seen. For small energies there is good agreement between 
the x-only approximation and the calculation by Foster 
et al. (while this agreement can already be seen in the 
one-fold ionization, it is not conclusive there because of 
the large uncertainties induced by the observable prob- 
lem). This is another hint (additional to the fact that 
/ < 1) that the x-only approximation is better suited for 
our system than the LK approximation. The reasons for 
the deviations between the results of our x-only calcu- 
lation and that by Foster et al. calculation for medium 
and high energies are not clear, but the good agreement 
of our results with the other shown data suggests that 
our results are solid in that region. 



While no reasonable bounds can be given for the two- 
fold ionization probabilities, the approximations (14), 
and (19) in (16) can be applied, yielding the results dis- 
played in Fig. 7. Once again, x-only and LK approxi- 
mations agree at high energies, while the TCS obtained 
from the LK approximation is higher at smaller energies. 
This confirms what has already been learned from the 
one-fold ionization. It can be seen that for two-fold ion- 
ization the observable problem is even more important 
than the question of the right approximation for the cor- 
relation potential. These large deviations introduced by 
the observable problem have to be considered a setback 
for the TDDFT approach. In order for TDDFT to be a 
viable approach for the calculation of two-fold ionization 
it would be necessary to find a much better approxima- 
tion for the probability as a functional of the density than 
that which we applied. 



So far, we investigated impact energies larger than 
5 keV. However, lower energies are in principle accessible 
experimentally (at CERN) and might prove to be inter- 
esting. Straightforward application of our method leads 
to a numerical problem concerning the scattering states: 
While they are constructed in a way that is suited to de- 
scribe whether an electron is in a scattering state or not, 
the dynamics of an electron in a scattering state are not 
described well. Especially, the average electron-target 
distance is strongly underestimated and the electrons are 
found to be in close vicinity of the target. When the den- 
sity obtained from such a state is used to calculate the 
Kohn-Sham potential, those electrons screen the target 
potential very strongly, leading to a much too high ion- 
ization probability. At sufficiently high energies this nu- 
merical effect is of no great concern since the time devel- 
opment is stopped before the effect grows large enough. 

For the results presented in this section we circumvent 
this problem by calculating the Kohn-Sham potential not 
from the full density, but only from that part obtained 
from the bound states: Let {|^}} be the bound basis 
states, {|^{)} the scattering states. The full density is 



n(r 1: t) 



(23) 



We define the bound density by 



£«i(*)Vi(ri 



(24) 



In the same way a free density can be defined (note, 
however, that in general n ^ + n b ). Calculating the 
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potential by using only n b is equivalent to the assumption 
that an electron in a scattering state moves away very 
quickly and has no further influence on the problem. 

In the low energy range, the straight-line approxima- 
tion for the projectile is no longer justified. Instead we 
use a classical trajectory determined by the initial kinetic 
energy and the approximate force by the target nucleus 
and the electrons in the following way: We assume the 
target (including the electrons) to be a point charge with 
the time-dependent effective charge 



Qeff(t)=2 



d 3 rm(n,0), 



(25) 



S(R(t)) 



where R(t) is the distance between the projectile and the 
target nucleus at time t. That means that the nuclear 
charge of 2 is reduced by the electrons that were initially 
closer to the target nucleus than the projectile is at the 
considered time [24] . It is not taken into account that the 
target system is no longer inertial when the force between 
target and projectile is nonzero. This would lead to an 
additional term in the TDSE [22]. In a classical picture 
this effect can be described as the nucleus being pulled 
away from under the electrons which remain in place due 
to their inertia. We therefore expect that consideration 
of this effect would magnify the ionization cross sections. 

For one-fold ionization (Fig. 8) as well as for two-fold- 
ionization (Fig. 9) we find that the curved trajectory 
leads to a higher ionization cross section. This is due to 
the fact that the effective impact parameter is reduced 
and the effective interaction time enhanced. While the 
quantitative results should be taken with a grain of salt 
due the simple modeling of the force on the projectile, 
it seems clear qualitatively that the trajectory effect is 
significant. We reiterate that non-inertial effects have 
been neglected, so that the real effect of the (classical) 
projectile-target interaction might be even larger. 
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FIG. 8: (Color online) Total cross section for one- fold ion- 
ization for small energies. Present results within x-only and 
binomial approximations and with straight-line and curved 
trajectories for the projectile motion. Experiment: Knudsen 
[13]. 



0.1 



co 
o 



0.01 



0.001 



x-only, binomial, curved 
x-only, binomial, straight 
Knudsen 



0.1 



1 

energy [keV] 



10 



FIG. 9: (Color online) Total cross section for two- fold ion- 
ization for small energies. Present results within x-only and 
binomial approximations and with straight-line and curved 
trajectories for the projectile motion. Experiment: Knudsen 
[21]. 



IV. CONCLUSIONS 



In this paper we have applied several approximations 
within TDDFT to the antiproton-helium collision sys- 
tem. We conclude that the correlation potential from 
Ref. [1] is not suited for this system because its main 
characteristic, the discontinuity at / = 1 does not come 
into play. Hence, until a better approximation is found, 
the x-only approximation should be preferred for calcu- 
lations within TDDFT. We saw that functional correla- 
tion and dynamic correlation are about equally large for 
one-fold ionization, whereas for two-fold ionization the 
functional correlation is the much more important effect. 
In the case of single ionization all those effects are neg- 
ligible for sufficiently high energies. In order to test dif- 
ferent approximations for the correlation potential, one 
should therefore focus on energies below 100 keV. Since 
functional correlation is also important in that regime, it 
is the average ionization that should be used for bench- 
marking. For energies below 5 keV we showed that there 
is a trajectory effect that significantly magnifies one- fold 
as well as two-fold ionization. 
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